11436 - Cubes EXTREME!/j.3.cpp

   1 #include<iostream>
   2 #include<math.h>
   3 #include <set>
   4
   5 using namespace std;
   6
   7 typedef unsigned long long ull;
   8
   9 inline ull cube(ull x){ return x*x*x; }
  10
  11 int main(){
  12   ull n;
  13   while(cin>>n &&n){
  14     bool found = false;
  15     for (ull k = 1; k*k<=n && k <= 29241; ++k){
  16       if (n % k != 0 ||
  17           3*k*k*k*k > 12*k*n){
  18         continue;
  19       }
  20       ull y = ( -3*k*k + (ull)sqrt(12*k*n - 3*k*k*k*k) ) / (6*k);
  21       if (y > 0 && cube(y+k) - cube(y) == n){
  22         //      cout << "k: " << k << endl;
  23         cout << y + k << " " << y << endl;
  24         found = true;
  25         break;
  26       }
  27       /*if (k*k < n){
  28         k = n/k;
  29         ull y = ( -3*k*k + (ull)sqrt(12*k*n - 3*k*k*k*k) ) / (6*k);
  30         if (y > 0 && cube(y+k) - cube(y) == n){
  31           //    cout << "k: " << k << endl;
  32           cout << y + k << " " << y << endl;
  33           found = true;
  34           break;
  35         }
  36         k = n/k;
  37         }*/
  38     }
  39
  40     if (!found){
  41       cout << "No solution" << endl;
  42     }
  43   }
  44   return 0;
  45 }